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Abstract. We investigate the formation of methane line at 2.3 
/im in Brown Dwarf Gliese 229B. Two sets of model parame- 
ters with (a) T cS = 940K and log(ff) = 5.0, (b) T cff = 1030K 
and log(g) = 5.5 are adopted both of which provide excel- 
lent fit for the synthetic continuum spectra with the observed 
flux at a wide range of wavelengths. In the absence of observa- 
tional data for individual molecular lines, we set the additional 
parameters that are needed in order to model the individual 
lines by fitting the calculated flux with the observed flux at 
the continuum. A significant difference in the amount of flux 
at the core of the line is found with the two different models 
although the flux at the continuum remains the same. Hence, 
we show that if spectroscopic observation at 2.3pm with a res- 
olution as high as R ~ 200, 000 is possible then a much better 
constraint on the surface gravity and on the metallicity of the 
object could be obtained by fitting the theoretical model of 
individual molecular line with the observed data. 
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broad wings of the K I and the Na I doublets and not due to 
the presence of condensates. All these model spectra together 
with the bolometric luminosity of the object (Leggett et al. 
1999) and the evolutionary sequences (Saumon et al. 1996) can 
constrain the effective temperature of the object very tightly, 
however the surface gravity is still poorly constrained. It is 
shown (Saumon et al. 2000) that a multi-parameter fit of the 
observed spectrum both for the K-band and for the red end is 
possible with different metallicities. Therefore, in order to de- 
termine the physical properties of the atmosphere of Gl 229B 
and the other T dwarfs uniquely, more comprehensive theoret- 
ical modeling of the observed spectrum is required. One of the 
most important studies is the individual line formation by the 
most abundant molecules, eg. CH4, H2O, NH3 etc. 

In this paper, we for the first time, attempt to model the 
line formed by methane at 2.3 /J,m and show that if individual 
lines of abundant molecules, in particular that of methane can 
be resolved observationally, then stringent constraint on the 
surface gravity, on the metallicity and on the temperature at 
the bottom of the atmosphere, where the optical depth is very 
high, can be obtained. 



1. Introduction 

The discovery of methane bands in the spectrum of Gl 229B 
(Nakajima et al 1995, Oppenheimer, Kulkarni, Matthews & 
Nakajima 1995, Geballe, Kulkarni, Woodward, & Sloan 1996) 
has not only helped to identify Gl 229B as a Brown Dwarf but 
also prompted the creation of the new spectral class of T dwarfs 
(Kirkpatrick et al 1999). The synthetic continuum spectra of 
the object has been obtained with and without dust particles 
(Marley et al. 1996, Griffith, Yelle & Marley 1998, for a review 
see Allard et al. 1997). Although the incorporation of conden- 
sates can explain the very rapid decline of the observed con- 
tinuum flux in the optical region, it is shown (Tsuji, Ohnaka, 
& Aoki 1999, Burrows, Marley, & Sharp 1999) that the pres- 
sure broadened red wing of the 0.77 fim K I doublet could 
also account for the observed features of the continuum flux 
shortward of 1.1 /im. The spectrum of the first field methane 
T dwarf SDSS 1624+0029 shows a broad band absorption fea- 
ture centered at 7700 A which is interpreted (Liebert et al. 
2000) as the K I 7665/7699 resonance doublet. Hence, it is 
most likely that the shape of the red spectrum is due to the 
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2. The Line Transfer Equations 

In order to model the individual line of a particular molecule 
one needs to solve the Non-LTE line transfer equations. The 
two level atom line transfer equation in the plane parallel strat- 
ified medium can be written as (Mihalas 1978) 

± / J( ^' Z) = Ku, z)[S(u, ± M , z) - I{v, ± M , *)] , (1) 

where jj, is the cosine of the angle made by the ray to the 
normal. The total absorption coefficient k(v, z) — ki (z)(j>(v, z) + 
k c (v,z), where ki is the frequency integrated line absorption 
coefficient, /3(z) = k c /ki, k c is the absorption coefficient for the 
continuum and ^(v, z) is the line profile function. The source 
function is given by 



S{v,±n,z) = 



{u,z)Si{u,±n,z)+P{z)Sc{z) 
<t>{v,z)+l3{z) 



(2) 



where S c is the continuum source function which is set 
equal to the Planck function B(v, z) and the line source func- 
tion Si is written as 



Si(u,±fi,z) = 



1 



(y ', z)dv x 



2 



R(v, n : u' , fi')I(p' , fjf , z)d// + eB(^, 2) . 



(3) 



where e is the probability per scattering that the photon is 
destroyed by collisional de-excitation and is given by (Mihalas 
1978) e ee e'/(l + e'); e' ee C* u( (1 - e hv/kT )/A ul , C ui is the rate 
of collisional de-excitation and A u i is the rate of spontaneous 
emission. R is the redistribution function defined according 
to a hybrid model prescription by Rees & Saliba (1982) and 
is given by R(v, fi; v' , //) = P(fi, fj,')R(v,u'), where P(/n, //) is 
the phase function and R(v, v') is the angle- averaged frequency 
redistribution function for isotropic scattering. 

3. The Absorption Coefficient 

Under chemical equilibrium among different atomic and molec- 
ular species in their standard states, the number density of 
CH4 for the molecular equilibrium reaction C + 2H2 CH4 
can be written as (Sharp 1985) 



ncH 4 — ncn H: 



( fr 2 A jpP e Qch 4 
\2-KmkT ) QcQ 2 



(4) 



H 2 



where Qc and Qh 2 are the partition functions for C and 
H2 , m is the "multiple" reduced mass and D is the disso- 
ciation energy of CH4 in eV, nc and n H2 are the number 
densities of C and H 2 , 9(z) = 5040/T(z) where T(z) is in 
Kelvin, and the other symbols have their usual meaning. Al- 
though chemical equilibrium of C and H2 should be coupled 
with some other molecules, dominance of CH4 molecule in the 
atmosphere would make the above equilibrium the most proba- 
ble one. Hence, we ignore, for the present, chemical equilibrium 
of C and H2 with other molecules. For the purpose of calcu- 
lating the number densities of the molecules in various energy 
levels, we assume that the lower atmosphere of the object is 
in local thermodynamic equilibrium (LTE). Then the number 
of particles nj, in a specified rotation level J is related to the 
total number of particles in all levels, ncH 4 , according to the 
following relation (Larson 1994): 



nj _ (2J+l)e- E -'/ feT M 



ncH 4 



Qch 4 



(5) 



where (2J+1) is the degeneracy factor and e~ Ej ^ kT is the 
Boltzmann factor. Now the integrated line absorption coeffi- 
cient can be written as 



h{z) = 



2 2 

fnj(z) = / 

m e c m e c 

,-Ej/kT(z) 



h 2 



2ixmkT{z) 



10 



D o 0(z) 



(2J+l)e- 



QcQ 2 



-n c (z)n H2 (z) . 



(6) 



H-2 



where / is the oscillator strength of the transition. We calcu- 
late the Q's for H2 and C using polynomial expressions given 
by Sauval & Tatum(1984). For a small value of J, the Boltz- 
mann factor reduces to 1. The dissociation energy D of CH4 
is taken to be 4.38 eV (Jorgenson 1994). We use the continuum 
absorption coefficient k c provided by D. Saumon. 



is clear and dominant. The monochromatic radiative transfer 
equations are solved numerically by using discrete space theory 
(Peraiah and Grant 1973). We have adopted two sets of model 
parameters, model (a): T o ff=940 K, log g=5.0 (g in cms -2 ) 
and [M/H]=-0.3 (K band) and model (b): T cff = 1030 K, log 
g=5.5 and [M/H]=— 0.1 (K band). The temperature and pres- 
sure profiles for both the models are obtained from M. Marley 
(private communication). The values of T e g and the surface 
gravity g are constrained by the bolometric luminosity of Gl 
229B and the evolutionary sequence of Saumon et al (1996). 
These synthetic spectra fit the entire observed spectrum of 
Gl 229B except in the near infra-red region. It should be men- 
tioned that a good fit with the observational data for the entire 
wavelength region is not possible with the same value of the 
metallicity and the surface gravity. The above sets of model pa- 
rameters are two of many optimal sets of parameters that can 
produce good fit. In order to match the observed flux shortward 
of 1.1 fim that declines very rapidly, either dust particulates 
or alkali metals have to be incorporated. However, the size of 
the grain that is needed to explain the observed spectrum in 
the optical region is too small to play any role in the infra- 
red region we are interested in the present work. Hence in this 
spectral region the law of Mie scattering that describes the an- 
gular distribution of photons in a dusty atmosphere reduces 
to that of Rayleigh scattering. Therefore, the synthetic spec- 
trum in this region matches well with the observed flux when 
Rayleigh's law of scattering is used. Recent ovservational evi- 
dence (Liebert et al. 2000) implies that there is no compelling 
reason to introduce dust or additional opacity source in the 
atmosphere of methane dwarf. Therefore in the present work 
we have not incorporated dust opacity. 

Figure 1 shows that there is no difference in the calculated 
continuum flux with the two sets of model parameters. Both 
the models fit the observed data very well. 
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Fig. 1. Synthetic continuum spectrum of Gl 229B: broken line is for 
the model (a) with T c ff = 940K and log(g) = 5.0 (g in cms -2 ), solid 
line for the model (b) with T eS = 1030K and log(g) = 5.5 



4. The Synthetic Continuum Spectra in the Infra-red 
Region 

We present the synthetic continuum spectra of Gl 229B (Fig- 
ure 1) in the infra-red region where the signature of methane 



5. Results and Discussion 

The continuum opacity and the temperature-pressure profile of 
the atmosphere are set by matching the synthetic continuum 
spectrum with the observed spectrum for the entire wavelength 
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region. The modeling of individual lines needs the values of ad- 
ditional parameters. In the absence of observational data for 
individual lines, we set the value of these parameters by match- 
ing the calculated flux in the continuum with the observed con- 
tinuum flux at 2.3 /jm. Since the line is very weak in intensity, 
we choose the profile function as (Mihalas 1978) 



-[(»-» )/A» D (^)] 2 



(7) 



where Aud is the thermal Doppler width and vq is the 
line center. It should be worth mentioning that individual 
molecular lines are usually not saturated enough so that pres- 
sure broadening is less important for them as compared to 
strong atomic lines. Moreover, molecular lines often overlap so 
strongly that their wings are completely masked (Schweitzer 
et al. 1996) and only the Gaussian line cores of the strongest 
molecular transitions are observed. The atmosphere of a brown 
dwarf is therefore only weakly sensitive to the Van der Walls 
damping constant. Nothing is known, at present, about the 
rotation of the brown dwarf Gl 229B around its own axis of ro- 
tation. If the projected velocity 'usini' of the object is greater 
than 2 to 5 kms -1 then rotational broadening could be signif- 
icant. However, in the present work we have neglected rota- 
tional broadening in order to make the results consistent with 
the calculation of the evolutionary sequences by Saumon et 
al. (1996) that constrains the surface gravity and the effective 
temperature of the object. The whole purpose of the present 
work is to show that with different values of the surface grav- 
ity and the metallicity, the flux at the line core is significantly 
different although it is the same in the continuum. Since rota- 
tional broadening would affect the spectrum equally for both 
the models, it is not important in the context of the present 
work. We assume complete frequency redistribution and use 
Rayleigh phase function for the angular redistribution. 

The parameters that are to be set in order to model the 
CHi line at 2.3 /an are e, /, nc{z), nH 2 ( z ), an d the degen- 
eracy factor (2 J + 1). We define s = (2 J + l)/nc(z)n| 2 («) 



guided by equation (|6j) and set the values of s and e such that 
the calculated value of the flux at the continuum matches with 
the observed continuum flux. We assume that e is indepen- 
dent of the geometrical depth. This is valid if the CH4 line 
formation is confined to a narrow region in the atmosphere. 
After testing several empirical laws we adopt the usual in- 
verse square law with respect to the geometrical depth for the 
variation of the number density of C and Hi. This should be 
verified when observational data becomes available. Experi- 
mental determination of the oscillator strengths for different 
transitions at 2.3 /im and at the relevant temperature could 
further constrain the number density of C and H2 and hence 
the abundance of methane in the atmosphere of Gl 229B pro- 
vided chemical equilibrium exists between the molecule and 
the atoms. We have solved the radiative line transfer equa- 
tions by using discrete space theory (Peraiah 1980). The the- 
oretical models for the methane line at 2.3 /jm are presented 
in Figure 2 . For the model (a) we find that the flux at the 
continuum matches with its observed value when e = 0.08 and 
s — (2J + l)fn\n% = 5.9 x 10 37 where ri\ and n,2 are the num- 
ber densities of C and H2 respectively at the bottom of the 
atmosphere. The moderately high value of e is consistent with 
the temperature at the lower atmosphere where the lines are 
formed. For the model (b) the flux at the continuum matches 
with the observed flux when e = 0.08 and s = 2.7 x 10 39 . Since 



the temperature of the lower atmosphere for the two models 
does not differ much, the value of e remains the same for both 
the models. However s differs substantially for the two models. 
This is because of the fact that for the model (a) the contin- 
uum opacity is less and therefore one has to reduce the line 
opacity in order to keep the right ratio (/3) between them so 
that the calculated flux matches with the observed flux at the 
continuum. However, the higher value of s in the model (b) 
makes the line opacity higher than that for the model (a) and 
so substantial decrease in the calculated flux at the line core 
is obtained. The difference in the flux reduces as we go to- 
wards the wings. It is worth mentioning that the wings could 
be masked by other lines whereas the continuum is overlapped 
by several lines. 




AX(A) 

Fig. 2. Emergent flux against wavelength from the line center 
(2.3/im): the curve 'a' is for the model (a) and the curve 'b' is for 
the model (b) 



Figure 2 shows that a spectral resolution as high as 200,000 
at 2.3/im is needed in order to investigate the individual molec- 
ular lines. This may be possible with an appropriate combi- 
nation of the telescope and the instrument. For example, the 
Cooled Grating Spectrometer 4 (CGS4) available in UKIRT 
(United Kingdome Infra-Red Telescope) has a spectral reso- 
lution upto 40,000. If observation of Gl 229B is possible at 
present by UKIRT with the maximum resolution power of 
CGS4 then keeping the signal to noise ratio (which is propor- 
tional to the diameter of the telescope and to the square root of 
the integration time of exposure) unaltered, a 10 m telescope 
(such as Keck I) can obtain the desired resolution by using a 
similar type of spectrometer provided the resolution of the in- 
strument is increased by about five times and the integration 
time of exposure is increased by about 2.5 times. 

It is found that the numerical values of e and s are very 
much sensitive to the emergent flux. The difference in the flux 
at the line core is clearly due to the different values of the 
surface gravity and the metallicity. Theoretical modeling of 
the continuum flux provides a few possible combinations of 
the metallicity, surface gravity and the effective temperature 
that are appropriate in explaining the observed continuum flux. 
The observational fit of the flux at the line core would decide 
which one of these combinations should describe the physical 
properties of the atmosphere. The physical parameters that are 
needed to model the individual molecular line will also be fixed 
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once observational data is available. Therefore, in conclusion 
we would like to emphasize that a theoretical fitting of the 
observed flux for the individual lines of any of the dominant 
molecules, especially methane, would not only provide a much 
better understanding on the abundance of that molecule and 
the temperature of the lower atmosphere but also improve the 
constraints on the value of the surface gravity of brown dwarfs. 
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